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Abstract 



The transition from order to chaos in atomic nuclei has been studied 
analytically and numerically using a quadrupole-quadrupole residual inter- 
action. This interaction leads to chaotic behaviour, but the critical energy 
Ec ^ 12.6 MeV, corresponding to the onset of chaos, is higher than that of 
the experimental one. 
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1. Introduction 

As is well known, in many-body systems states of different " complex- 
ity]^' generally coexist [1,2,3]. In particular, in atomic nuclei the "hierarchy" 
of complexity can be represented by a simple picture. In zero order ap- 
proximation the relevant elementary excitations (single-particle, rotations, 
vibrations) may regarded as independent modes and then the interaction 
between these elementary modes can be considered [4]. At the end of the 
chain of complexity there are the so called "stochastic" or "chaotic states", 
whereas a single mode can be considered "regular". 

The experimental data of nuclear spectroscopy (see, e.g. references [5,6]) 
suggest regular states near the ground state (0-3 MeV) and chaotic states 
near the neutron emission threshold (6-8 MeV). On the other hand, in the 
framework of the so called "rotational damping", the authors of references 
[7] have shown that the surface delta interaction produces chaotic behaviour 
while a pairing plus quadrupole force is unable to do so. The aim of this 
paper is to analyze in detail the quadrupole-quadrupole residual interaction 
showing that it is able to reproduce the transition order-chaos in the atomic 
nuclei, but at an energy (~ 12.6 MeV) higher that that of the experimental 
one (~ 8 MeV). 

The structure of the article is as follows: in section 2 we present a 
schematic nuclear model that describes two nucleons in a mean-field in- 
teracting with a separable quadrupole-quadrupole residual interaction; in 
^In this paper the terms "complexity" and "hierarchy" have been used in a general 
sense; for precise definitions, sec [17]. 
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section 3 we apply the curvature of potential energy criterion to calculate 
the energy Eq of the onset of chaos and in section 4 we plot the Poincare 
sections of the model for bound energies. Finally in section 5 we diagonalize 
the quantum hamiltonian and calculate the distribution P{S) of spacings 
between adjacent levels. 

2. The schematic quadrupole-quadrupole interaction 

In our model nucleons move independently in an oscillator potential under 
the influence of a separable quadrupole-quadrupole effective residual inter- 
action. The hamiltonian of the system is given by [8]: 

i=i i=i 



ie T/ A 2 

E {-^r{rlr.MA^)}{r]Y,_^{e,A,)}. (1) 



• 

i<j m=—2 



where Rq — 1/(3 is the range of the interaction and A the number of nucleons. 
We consider only 2 nucleons and suppose that they move in one dimension, 



so that = Y (^)^m,o, and we have: 

H = ^(pl + pI) + V,P\l + rl) - (2) 
To simplify the problem we perform the transformation: 



/ 1 /ma; 

where ma;^/2 — Vo/i?o- So the hamiltonian becomes: 

H = e{pi +P2 + qi + q2) - Xolfll, (4) 



with € = hu) and x = S/3{h'^ /mRl). 

In nuclei e = /lo; ~ 41 A~^/^ MeV is the single particle energy for the 
harmonic oscillator potential [8] and Rq — 1/P 1.2A^/^ fm is the nuclear 
radius, and so we have x — 74A"^/^ MeV. 

3. The onset of chaos with the curvature criterion 

We apply the criterion of curvature [9] to the hamiltonian (4) in order to 
calculate the critical energy of the transition to chaos, which will be indicated 
by Ec- The criterion of curvature is based on the estimation of the rate of 
separation of neighbouring trajectories in the phase space. To calculate the 
time evolution of the dynamical system with hamiltonian: 

H ^ ^{pi+P2) + V{n,r2) (5) 
the following equations have to be solved: 

dt^ dp ' dt'^ df ' 
where r = (ri,r2) and p = {pi,P2)- The deviations from the two initially 
neighbouring trajectories {f{t),p{t)) are given by: 

df{t) =m- r{t) , 6p{t) =m- Pit) (7) 
and the linearized equations of motion for the deviations are: 

= M-'dp{t), = -S{t)dr{t) (8) 

where M^^ = Sijm~^, and: 



The stability of the dynamical system is then determined by the eigenvalues 
of the 4x4 matrix: 

/ \ 
m={ ■ (10) 

\-s{t) ; 

If at least one of the eigenvalues Xi{t) is real, then the separation of the trajec- 
tories grows exponentially, and the motion is unstable. Imaginary eigenvalues 
correspond to stable motion. To diagonalize the matrix r{t), we must first 
solve the equations of motion of the differences (8). The problem can be 
significantly simplified by assuming that the time dependence can be elimi- 
nated, i.e. r{f{t)) — r(f). The eigenvalues then are: 

Al,2,3,4 = ±M±\/&^^]^ (11) 

where: 



drl dri dridr2 

Now, if 6 > then with c > the eigenvalues are purely imaginary and 
the motion is stable, meanwhile with c < the pair of eigenvalues become 
real, and this leads to exponential separation of neighbouring trajectories, 
i.e. chaotic motion. The parameter c has the same sign as the Gaussian 
curvature K{ri, r2) of the potential-energy surface: 

dr'f 



Let us now return to our nuclear problem. The potential energy is: 

V{ri,r2) = e{rf + rl) - xrjrl; (15) 



it has one minimum for ri = r2 = with = MeV, and four saddle points 

for n = ±yf , T2 = ±yj, and n = i^J, r2 = T^J, with = f 23.3 

MeV. These are the points for which the curvature criterion is exact (Vi? = 

0); the origin is a stable elliptic point with Ai 2,3,4 = ±y^2)i and the four 

saddle points arc unstable hyperbolic points with Ai,2 = ±2 and A3,4 = ±2i. 

The equipotentials (curves of constant potential V) are shown in figure 

-, ,/-) = — and 
x' V X ^ 
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hounded motion ior E < ^. We concentrate here upon the region < £■ < ^. 

At low positive energies, the motion near the minimum of the potential 
energy, where the curvature is positive, is periodic or quasi-periodic and is 
separated from the region of instability by a line of zero curvature. If the en- 
ergy is increased, the system will, for certain initial conditions, be in a region 
of negative curvature where the motion is chaotic. In accordance with this 
scenario, the energy of order^chaos transition Ec is equal to the minimum 
value of the line of zero gaussian curvature K{ri,r2) of the potential-energy 
surface of the system. For our potential, the gaussian curvature vanishes at 
the points that satisfy the equation: 

2 

4e(l - ^rl - ^rl - S^r^r^) = 0. (16) 



The energy on the zero-curvature line is determined by the expression: 

ViK = 0, r2) = - Jf7}_, + + erl (17) 



(see also fig. 2). It is easy to show that the minimal energy on the zero- 
curvature line is given by: 

K.,„(X = 0,r-i) = ^, (18) 



and occurs at r^i = ±-0- This is the critical energy of the transition to chaos 
Ec of the model, and we have Eq — 12.62 MeV. We observe that the energy 
of the onset of chaos {Ec ~ 12.6 MeV) is much higher than that required for 
the atomic nuclei one (chaotic behaviour ~ 8 MeV) [5,6]. 

In figure 3 we show the rate of regular points of configuration space ob- 
tained with the curvature criterion; the transition to chaos seems to be quite 
smooth. 

4. Classical calculations 

The curvature criterion is only able to characterize the local behaviour 
of the system (for example local instability) and gives a rough signature of 
the global properties (e.g. global instabihty) [10]. As is well known, global 
properties are more appropriately studied using Poincarc sections [12,13]. 

The classical equations of motion obtained from the hamiltonian (4) are: 

ri = 2epi, 7-2 = 2ep2 
Pi = -2eri + 2xrirl, p2 = -2er2 + 2xrlr2- (19) 

We used a fourth order Runge-Kutta method [11] to compute the classical 
trajectories. Poincare sections [12,13] have been produced for a variety of 
energies, showing the transition from almost totally regular motion at low 
energies to almost totally irregular motion at high energies. Conservation of 
energy restricts any trajectory in four-dimensional phase space to a three- 
dimensional energy shell. At a particular energy, therefore, the restriction 
ri = defines a two-dimensional surface in phase space. Each time a par- 
ticular trajectory passes through the surface, i.e. each time it crosses the r2 
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point is plotted at the position of intersection (r2,p2)- We employ a 
first-order interpolation process to reduce inaccuracies due to the use of a 
finite step length [13]. 

Regular regions on the surface of section plots are characterised by sets 
of invariant intersection points. The surface of section pictures (fig. 4) show 
that motion is almost wholly regular for E = 15 MeV. However, for E = 21 
MeV we see evidence of bifurcation: one can count at least 6 new elliptic 
points and some points are distributed irregularly. For E — 23 MeV the 
number of irregularly distributed points increases but many regular regions 
remain. 

These numerical results are in good agreement with the smooth transition 
to chaos predicted by the curvature criterion. 

5. Quantum calculations 

In quantum mechanics one cannot apply classical concepts and methods 
directly since the notion of trajectory is absent. Nevertheless, many efforts 
have been made to establish the features of quantum systems which reflect 
the qualitative difference in the behaviour of their classical counterparts [2]. 
Many schematic models [2,3] have shown that this difference reveals itself in 
the properties of fluctuations in eigenvalue sequences. The spectral statistics 
for the systems with underlying chaotic behaviour agree with the predictions 
of the random matrix theory. By contrast, quantum analogs of classically 
integrable systems display the characteristics of Poisson distribution. 

To obtain the quantum mechanical energy levels of the system we can 
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write the hamiltonian (4) with the creation and destruction operators: 

ttk = [r-k + ipk), al = {vk - ipk), (20) 

where k — 1,2, so we have: 

H = e{atai + of 03 + 1) - ^^(01 + 0^)^(02 + atf. (21) 

The eigenvalues can be calculated by diagonalizing the hamiltonian in the 
basis \n1n2 > of the occupation numbers of the two harmonic oscillators. 
The matrix elements can be written: 

< n[n'2\H\nin2 >= e(ni + ns + ^)Sn[ni^n2n2 + 

+Y^[Vni(ni - 1)<5„;„^_2 + ^ {n, + l){n, + 2)5^,^^^^, + {2n, + l)5„,„Jx 

X [^n,{n, - 1)S^>^^_, + V(ri2 + l)(n2 + 2)5„.„^^, + (2^2 + 1)5„.„J (22) 

and each submatrix can be labelled by the parity of the occupation numbers 
Til, n2- By performing the unfolding procedure described in detail in refer- 
ence [14] , each spectrum has been mapped into one with quasi-uniform level 
density. 

The distribution P{S) of spacings between adjacent levels has been cal- 
culated and compared to the Brody [15] distribution: 

P{S) = a{q + l)S'^exp{-aS'^+^), (23) 

with: 

a^{T[^]r+\ 0<q<l. (24) 
0+1 
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The distribution interpolates between the Poisson distribution {q = 0) of 
integrable systems and the Wigner-Dyson distribution {q — 1) of chaotic 
ones. In figure 5 (a) the distribution P{S) obtained with all the four classes 
is plotted for < < 12 MeV with g = 0, and in (b) for 12 < < 23 MeV 
with q = 0.15. 

6. Conclusions 

In this paper we have shown that the quadrupole-quadrupole residual 
interaction leads to the transition from order to chaos. However the energy 
of the onset of chaos {Eq ~ 12.6 MeV) is much higher than the experimen- 
tal one of atomic nuclei (chaotic behaviour ~ 8 MeV) [5,6]. Therefore the 
quadrupole-quadrupole force is unable to reproduce experimental nuclear 
data and higher multipole components should be added [16]. 
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Figure Captions 



Figure 1: The potential energy of the model {A = 200). 

Figure 2: The contour plot of the potential energy and the region of negative 

curvature (marked by •). 

Figure 3: Ratio of stable points in phase space obtained from the curvature 
criterion. 

Figure 4: The Poincare sections of the model: {&) E = 7 MeV, (b) E = 15 
MeV, (c) E = 23 MeV. 

Figure 5: The distribution P{S) of spacings S between adjacent levels: (a) 
< £; < 12 MeV (q = 0), (b) 12 < < 23 MeV (q=0.15). 
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